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I.  INTRODUCTION 


Large  space  structures  have  been  identified  as  having  a  number  of 
potential  applications  for  national  security  and  future  energy  needs.  Typical 
applications  include  solar  power  station,  large  space  mirror,  large  space 
antenna,  space-based  radar,  and  multipurpose  large  space  platform.  Various 
missions  which  require  large  lightweight  structures  and  the  major  design 
requirements  of  such  structures  have  been  summarized  by  NASA  in  Refs.  1  and  2. 
Although  prioritized  specific  missions  are  yet  to  be  defined,  several  novel 
and  innovative  design  concepts  have  been  suggested  for  future  applications 
(Refs.  3  through  14).  These  generic  conceptual  ideas  provide  valuable  informa¬ 
tion  regarding  the  technology  needed  for  developing  structurally  efficient  low 
cost  systems.  The  size,  design  environment,  manufacturing  methods,  and  other 
characteristics  associated  with  large  space  structures  dictate  that  the  first 
time  such  a  structural  system  assumes  its  operational  configuration,  it  will 
occur  in  orbit,  thus  creating  unique  design  and  testing  problems.  Furthermore, 
the  cost  involved  in  deploying  such  structural  assemblies  to  their  full  capa¬ 
city  requires  a  high  level  of  confidence  in  analytical  methods  and  modeling 


l"Outlook  for  Space,"  NASA  Task  Group,  NASA  SP-386,  January  1976. 

^Hedgepeth,  John  M.,  "Survey  of  Future  Requirements  for  Large  Space  Struc¬ 
tures,"  NASA  CR-2621 ,  1975. 

^Woods,  A.  A.,  Jr.,  "Offset  Wrap  Rib  Concept  and  Development  (LMSC),"  Large 
Space  Systems  Technology,  NASA  Conference  Publication  2118,  1979. 

^Archer,  J.  S.,  "Advanced  Sunflower  Antenna  Concept  Development  (TRW),"  Large 
Space  Systems  Technology,  NASA  Conference  Publication  2118,  1979. 

^Montgomery,  D.C.  and  L.D.  Skides,  "Development  of  Maypole  (Hoop/Column) 
Deployable  Reflector  Concept  for  Large  Space  Systems  Applications  (Harris 
Corp.),"  Large  Space  Systems  Technology,  NASA  Conference  Publication  2118, 
1979. 

^"Modular  Reflector  Concepts  Study  (GDC),"  Large  Space  Systems  Technology, 
NASA  Conference  Publication  2118,  1979. 
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techniques  for  predicting  their  structural  responses  to  the  natural  and  induced 
environments.  Although  significant  emphasis  has  been  placed  on  developing 
design  concepts,  construction  methodology,  and  deployment  techniques  for  large 
space  structures,  extremely  limited  efforts  have  been  devoted  tp  developing 
efficient  analytical  methods  for  structural  modeling.  As  currently  envisioned, 
a  strong  candidate  for  a  large  low-mass  and  high-stiffness  structure  is  an  open 
truss  configuration.  The  design,  analysis,  and  testing  of  such  a  structure, 
with  a  large  number  of  structural  elements,  present  several  challenging  prob¬ 
lems.  Some  such  problems  are  to  develop  capabilities  for  conveniently  assess¬ 
ing  preliminary  designs,  efficiently  carrying  out  tradeoff  studies,  confidently 
predicting  structural  responses,  suitably  designing  ground  test  methodology, 
and  judiciously  recommending  space  verification  test  programs. 


7"DOD/STS  On-Orbit  Assembly  Concept  Design  Study  (GDC),"  SAMSO-TR-78-128 , 

1978. 

8"D0D/STS  On-Orbit  Assembly  Concept  Design  Study  (MMC) Martin  Marietta 
Corporation,  Report  MCR-78-113,  1978. 

^Agan,  W.  E. ,  "Erectable/Deployable  Concepts  for  Large  Space  System  Techno¬ 
logy  (Vought  Corp.),"  Large  Space  System  Technology,  NASA  Conference  Publi¬ 
cation  21ld,  1979. 

^Britton,  W.  R.  and  J.  D.  Johnston,  "Space  Spider  -  A  Concept  for  Fabrica¬ 
tion  of  Large  Space  Structures,"  AIAA  Conference  on  Large  Space  Platforms: 
Future  Needs  and  Capabilities,  September  1978. 

^Stokes,  J.  W.  and  E.  C.  Pruett,  "Structural  Assembly  in  Space,"  Large  Space 
Systems  Technology,  NASA  Conference  Publication  2118,  1979. 

12Nein,  M.  E.  and  F.  C.  Runge,  "Science  and  Applications  Space  Platforms," 

AIAA  81-0458,  AIAA  Second  Conference  on  Large  Space  Platforms,  February  1*81. 

^Johnson,  R.  R. ,  H.  Cohan,  and  G.  G.  Jacquemin,  "A  Concept  for  High  Speed 
Assembly  for  Erectable  Space  Platforms,"  AIAA  81-0446,  AIAA  Second  Confer¬ 
ence  on  Large  Space  Platforms,  February  1981. 

l^Stoll,  H.  W. ,  Systematic  Design  of  Deployable  Structures,"  AIAA  81-0444, 

AIAA  Second  Conference  on  Large  Space  Platforms,  February  1981. 
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cions,  which  is  boch  time-consuming  and  expensive.  The  discrete  field  approach 
taxes  advantage  of  the  regularity  of  the  structure.  It  is  an  extension  of  the  i 

classical  method  used  for  continuum  problems  in  which  equilibrium  and  compati¬ 
bility  equations  at  each  node  (or  joint)  of  the  structural  frame  are  formulated 
first  and  then  the  resulting  governing  equations  for  the  total  system  are 
solved  using  finite  differences  and  differential  calculus. 

The  state  of  the  art  of  the  discrete  field  method  is  summarized  in  Ref. 

15,  whereas  some  example  problems  illustrating  this  approach  are  given  in  Refs. 

16  through  18.  Although  the  discrete  field  method  is  useful  for  moderate-sized 
problems  with  simple  geometry,  it  becomes  increasingly  complex  when  applied  to 
large  intricate  configurations.  - 

| 

t 

t 


l^Dean,  D.  L.  and  R.  R.  Avent,  "State-of-the-Ar t  of  Discrete  Field  Analysis  ; 

of  Space  Structures,"  Proceedings  of  Second  International  Conference  on  f 

Space  Structures.  Edited  by  W.  J.  Supple,  University  of  Surrey,  Guildford,  j 

England,  September  1975. 

i 

10Renton,  J.  D. ,  "The  Related  Behavior  of  Plane  Grids,  Space  Grids,  and  | 

Places,"  Space  Structures,  Blackwell,  Oxford,  England,  1967.  j 

i 

l^Renton,  j.  D,  "General  Properties  of  Space  Grids,"  International  Journal  of 
Mechanical  Sciences,  Vol.  12,  1970. 

l^Dean,  D.  L.  and  C.  P.  Ugarte,  "Field  Solutions  for  Two-Dimensional  Frame¬ 
work,"  International  Journal  of  Mechanical  Sciences,  Vol.  10,  1968. 
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The  chird  approach  is  a  technique  to  replace  discrete  structural  frame¬ 
works  with  equivalent  continuum  models  (kefs,  iv  tnroucr.  2  5).  It  involve*  tr.e 
determination  of  equivalent  elastic  and  dynamic  properties  of  the  truss  struc¬ 
ture.  The  response  of  a  given  structure  is  then  predicted  by  solving  the  con¬ 
tinuum  model  under  similar  loading  and  boundary  conditions.  The  continuum 
problem  can,  in  general,  be  simplified  by  making  certain  logical  cinematic 
assumptions.  For  example,  a  large  plat f orm-tvpe  truss  could  be  modeled  as  an 
equivalent  continuum  plate,  whereas  a  long  truss  boom  could  be  represented  by 
an  equivalent  continuum  beam.  This  approach,  when  applied  to  large  space 
structures,  has  great  versatility  in  efficiently  determining  the  overall 
response  of  the  structure  to  the  induced  loading  as  well  as  in  carrying  out 
the  design  feasibility  and  tradeoff  studies.  Although  much  work  has  yet  to  be 
done  toward  accomplishing  "perfection,"  the  present  study  represents  a  step 
forward  in  the  state  of  the  arc  of  developing  modeling  methodology  for  large 
truss-type  structures  based  on  the  equivalent  continuum  approach. 


l^Mikulas,  m.  M. ,  Jr.,  H.  G.  Bush,  and  M.  F.  Card,  "Structural  Stiffness 
Strength  and  Dynamic  Characteristics  of  Large  Tetrahedral  Space  Truss 
Structures,"  NASA  TMX74001 ,  1977. 

^Bush,  H.  G. ,  M.  M.  Mikulas,  Jr.,  and  W.  L.  Heard,  "Some  Design  Considera¬ 
tions  for  Large  Space  Structures,"  Proceedings  of  the  AIAA/ASME  18th  Struc¬ 
tures,  Structural  Dynamics,  and  Materials  Conference,  Vol.  A,  1977. 

21Renton,  J.  D.,  "On  the  Gridwork  Analogy  for  Plates,"  Journal  of  Mechanics, 
Physics,  and  Solids.  Vol.  13,  1965. 

22Flower,  w.  R. ,  and  L.  C.  Schmidt,  "Analysis  of  Space  Truss  as  Equivalent 
Plate,"  Journal  of  the  Structural  Division,  Vol.  97,  ASCE,  1971. 

2^Heki,  K.,  "On  the  Effective  Rigidities  of  Lattice  Plates,"  Recent 
Researches  of  Structural  Mechanics,  Tokyo,  1968. 

2^Sun,  C..  T. ,  and  T.  Y.  Yang,  "A  Continuum  Approach  Toward  Dynamics  of 
Gridworks,"  Journal  of  Applied  Mechanics,  Vol.  91,  ASME,  1965. 

2^Nayfeh,  A.  H.  and  M.  S.  Hefzy,  "Continuum  Modeling  of  Three-Dimensional 
Truss-Like  Space  Structures-,"  AIAA  Journal,  Vol.  16,  No.  8,  1979. 


The  procedure  for  deriving  che  effective  elastic  and  dynamic  properties 
of  an  equivalent  continuum  modal  is  giver,  in  Section  TI,  uncle  the  analyte  to' 
solution  of  the  equivalent  continuum  mocel  is  discusses  ir.  Section  III.  .-.n 
example  problem,  for  verifying  accuracy  of  che  model  is  illustrated  in  Section. 
IV  followed  bv  the  conclusions  of  this  studv. 


II.  EQUIVALENT  CONTINUUM  MODEL 


A.  DETERMINATION  OF  EQUIVALENT  EFFECTIVE  PROPERTIES 


The  procedure  for  determining  the  effective  properties  of  an  equivalent 
continuum  model,  similar  to  that  of  Nayfeh  and  Hefzy  (Ref.  26),  was  indepen¬ 
dently  derived.  The  approach  is  straightforward.  For  a  given  truss  structure, 
all  sets  of  parallel  members  are  identified.  The  direction  cosines  for  each 
set  of  parallel  members  are  determined  with  respect  to  the  global  reference  co¬ 
ordinate  system.  Unidirectional  stiffness  properties  of  each  set  are  derived 
next  by  "smearing”  the  stiffness  of  individual  members  in  the  set  over  the 
effective  area.  Finally,  through  orthogonal  transformations,  contributions  of 
each  set  to  the  effective  continuum  properties  are  derived. 


For  a  linear  elastic  body,  the  stress-strain  relations  are  given  by 

*  1,2,3 


o.  . 

ij 


Cijkf  €ki 


(1) 


where  and  are  the  components  of  stress  and  strain  tensor,  respec¬ 
tively,  and  C^^j  are  the  elastic  constants.  Equation  (1),  in  matrix  form, 
is  written  as 
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26jjayfeh,  A.  H.  ,  and  M.  S.  Hefzy,  "Continuum  Modeling  of  the  Mechanical  and 
Thermal  Behavior  of  Discrete  Large  Structures,"  AIAA  80-0679,  AIAA/ASME/ASCE/ 
AHS  21st  Structures,  Structural  Dynamics,  Materials  Conference,  May  12-14, 
1980. 
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An  expression  for  determining  C.^^-  for  an  equivalent  continuum  in 

terms  of  the  geometric  and  material  properties  of  the  original  structure  shall 

now  be  derived.  If  (  8  .  .)  are  the  direction  cosines  of  the  mtn  parallel 

1  j  m 

member  set  and  (C1  )  are  the  smeared  stiffness  properties  of  that  set.  the 

pqrs  m 

global  contributions  (C..t  ,)  of  this  set  to  the  effective  equivalent  continuum 
°  ljki  m 

properties  are  given  by 


(c. ..  .)  =  (  8-  8-  jSu  8c  >  (c*  ) 

s  ljkx  m  ^ip  ^jq  ^kr  ns  n  pqrs 


ip  '"jq 


pqrs  m 


(3) 


where 


(8-  )  - 

ij  m  \  3Xj  J m 


U) 


or  alternatively 


(  8-  •)  =  [cos(x!  ,  x.)] 

'ij  m  i  j  m 


(5) 


where  x^  and  xf  are  the  global  and  local  coordinates,  respectively. 

Equation  (3)  can  be  further  simplified.  Because,  in  truss  structures, 
the  members  have  only  one  unidirectional  elastic  property  (the  modulus  of 

elasticity  along  the  principal  direction  of  the  truss  member)  Eq.  (3)  can  be 
rewritten  as 

(C...J  *  (  /S-,  £•,  Buy  £,,>  <6> 

ljki  m  ^il  ^jl  A'kl  ^fl  m  1111  m 


B. 


DETERMINATION  OF  (C '  ,  ) 

1111  m 


The  effective  unidirectional  elastic  property  of  a  set  of  parallel 

rods,  Fig.  1,  spaced  a  apart  is  given  by 

m 


(C*  .) 

1111  m 


E  A 
m  m 

a  h 
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(7) 
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and  Che  effective  mass  density  is  given  by 


a  h 

m 


where  A  is  the  cross-sectional  area,  E  is  the  modulus  of  elasticity,  and 
m  m 

is  the  mass  density  of  the  members  of  the  m  parallel  member  set.  The 
linear  dimension  h  is  a  reference  measurement  introduced  for  the  purpose  of 
obtaining  area  average. 


Fis.  1.  A  Set  of  Parallel  Rods 


The  equivalent  continuum  elastic  properties  are  now  obtained  by 

summing  the  contributions  of  all  n  parallel  member  sets 


C .  ...  *  (C.  .  ) 

ljkf  jr1.  ijkl  m 


and  the  effective  mass  density  is 


Z  »L 


C.  APPLICATIONS 


1.  Two-Dimensional  Framework 


The  procedure  for  obtaining  equivalent  continuum  properties  is  illus¬ 
trated  through  example  problems.  A  two-dimensional  framework  is  shown  in 
Fig.  2.  One  can  easily  identify  four  sets  of  parallel  members.  The  direction 
cosines  for  these  four  sets,  with  respect  to  x^,  Xj  axes,  are 
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Two-Dimensional  Framework 
(0°,  90° ,  +45°  Arrangement) 


If  A  ,  E  ,  and  p  (m  =  1  through  4)  are  the  cross-sectional  area, 
mm  m 

modulus  of  elasticity,  and  mass  density,  respect ively ,  of  the  members  in  these 

four  sets,  the  effective  mass  and  elastic  properties  of  each  set  are 
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and  all  other  (Cl  '  vanish. 
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When  Eqs.  (13)  are  substituted  into  Eqs.  (9),  the  equivalent  elastic 
continuum  properties  obtained  are 


f  h  tlAl  +  8  tE3A3  +  E4A45 


£  o 


C  ■  *-  (E  A  -  E  A  ) 

1122  8ih  U3  3  4A4; 


C1112  “  8fh  (E3A3  "  E4A4) 


7h  |E2A2  *4-  <E3A3  *  E4A4> 


C2212  ~  8fh  (E3A3  "  W 


and  the  effective  mass  density  is 


P0  "  fh  [  P1A1  *  P2A2  +  2  (P3A3  +  P4A4}] 


This  example  problem  demonstrates  the  simplicity  with  which  the  equiva¬ 
lent  continuum  effective  properties  can  be  derived.  The  method  shall  now  be 
applied  to  a  three-dimensional  truss. 
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A  Three-Dimensional  Problem  -  Tetrahedral  .russ 


A  sample  tetrahedral  truss  is  shown  ir.  Pis.  3.  It  consists  of  :vj 
parallel  layers  of  equilateral  triangular  meshes  connected  by  bracing  bars. 
Tne  two  triangular  paraiiei  layers  rorm  u  ,  +•  cO  arrangement  or  tee  center 
orientation.  The  geometric  and  material  properties  -f  the  truss  me-ders  are 

summarized  in  Table  1. 

If  h  is  the  height  of  Che  truss,  l  is  the  length  of  the  members  in  the 
top  and  bottom  layers,  and  d  is  the  length  of  bracing  bars  one  can  write. 


d 


(16) 


A  repeating  element  of  this  truss  is  shown  in  Fig.  4.  Six  sets  of 
parallel  members  are  identified  by  numbers  1  through  6.  The  direction  cosines 
for  these  six  sets  are  given  in  Table  2. 
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yj  *  /3d 


h/d 


Following  Che  example  problem  discussed  earlier,  c'ne  unidirectional 
effective  properties  for  the  members  in  the  top  and  bottom  layers  of  the 
tetrahedral  truss  can  be  determined  easily.  The  arranpemer.t  of  the  members  in 

O  O  -  . 

the  present  case  is  u  ,  *  cO  .  Table  o  lists  ere  efrective  unidirectional 
properties  of  the  three  sets  of  the  members  in  the  top  and  bottom  layers  of 
the  truss. 


Table  3.  Unidirectional  Properties  of  the  Members 
in  the  Top  and  Bottom  Layers 


Parallel  Member 
Set  Number 


Effective  Unidirectional  Properties 


Top  Layer 


Bottom  Layer 


Elastic 


Elastic 


1 ,  2 ,  ana  3 


yr  eiai 


2  s/T  Plkl 


2  JJ  E2A2  2  yT  P2A2 

3  fh  3  f  h 


The  effective  unidirectional  properties  for  the  bracing  members  are 
determined  next.  A  projection  of  the  tetrahedral  truss  on  a  plane  perpen¬ 
dicular  to  the  parallel  member  Set  4  (see  Fig.  4)  is  shown  in  Fig.  5.  In  this 
figure,  bold  dots  represent  the  members  of  Set  4,  whereas  the  effective  area 
occupied  by  each  member  is  shaded.  Unidirectional  effective  elastic  and  mass 
properties  for  the  Set  4  members  are  now  computed  by  smearing  the  stiffness  of 
each  member  over  the  effective  area.  Thus 


effective  stiffness  * 


~2  Ed  Ad 


effective  mass 


1  Pd  Ad 
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Equations  1.6),  (9),  and  (10)  are  now  used  for  deriving  the  elastic  and 
mass  properties  of  an  equivalent  continuum  model.  For  the  purpose  of  illus¬ 
tration,  expression  for  the  elastic  constant  is  derived  in  detail. 

From  Eq.  (6),  using  Tables  2  and  3  and  Eq.  (17) 
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Equation  (9)  yields 


Similarly,  other  C...  ,'s  can  be  derived: 
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and  che  equivalent  mass  density  is  given  by 


2  y/3 
fh 
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P2A2 


I  Vd 


( 19) 


This  completes  the  methodology  for  obtaining  the  equivalent  elastic  and  mass 
properties  for  the  continuum  in  terms  of  the  geometric  and  material  properties 
of  a  given  tetrahedral  truss. 

In  order  to  solve  the  continuum  problem  for  predicting  the  structural 
response  of  the  original  truss,  the  governing  field  equations  are  (Ref.  27): 

Equations  of  Motion 


F. 

l 


Constitutive  Equations 


Cijkf 


i.j.k.f  =  1,2,3 


(20) 


(21) 


^Fung,  Y.C.,  Foundations  of  Solid  Mechanics,  Prentice-Hall,  Inc.,  Englewood 
Cliffs,  NJ,  1965. 


31 


Strain-Displacement  Relations 


(.  . 

J-J 


(22) 


where  are  the  body  forces  and  u^  are  the  displacement  components.  The 
problem  in  this  form  is  so  wrought  with  mathematical  complexities  that  a 
general  solution  is  almost  impossible.  However,  it  is  often  possible  to  make 
the  problem  tractable  and  obtain  a  reasonable  solution  through  certain 
judicious  assumptions. 
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III.  SOLUTION  OF  THE  CONTINUUM  MODEL 


The  mathematical  difficulties  involved  in  determining  solutions  for  the 
equivalent  continuum  model,  derived  in  Section  II,  may  be  circumvented  by 
making  some  simplified  kinematic  assumptions.  For  example,  in  the  case  of  a 
tetrahedral  truss  representing  a  hypothetical  large  space  platform,  its  equi¬ 
valent  continuum  model  can  be  viewed  as  a  plate  (Ref.  28).  This  is  precisely 
the  approach  adopted  in  the  present  study. 

Since  the  displacement  components  vary  linearly  along  the  members  of 
the  truss,  it  is  assumed  that  the  displacement  components  u^,  u^.  and  u^  of 
the  equivalent  continuum  plate  model  vary  linearly  with  the  coordinate, 
i.e. 


u2  *  ‘u2  *  x3  *2 

0  0 
u3  *  u3  +  x3€33 


(23) 


where,  as  shown  in  Fig.  4.  n  ,  u_  ,  and  u,  are  the  mid-surface  displacement 

±  s  0 

components  (at  x^  *  0),  and  are  the  rotation  components,  and  is 

the  transverse  normal  strain.  The  mid-surface  displacements,  the  rotation 

components,  and  the  transverse  normal  strain  are  all  assumed  independent  of 

the  x^  coordinate. 


2®Szilard,  R. ,  Theory  and  Analysis  of  Plates,  Prentice-Hall,  Inc.,  Englewood 
Cliffs,  NJ,  1974. 


As  a  consequence  of  chis  displacement  field  assumption,  the  strain 
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(24) 


1 


The  stress  and  moment  resultants  for  the  equivalent  continuum  plate  are 
defined  as  follows: 


Stress  Resultants 
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Moment  Resultants 
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The  equivalent  stiffnesses  (viz.,  extensional,  coupling,  and  bending), 
for  the  continuum  model  are  as  follows: 

Extensional  Stiffness 
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Coupling  Stiffness 
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Bending  Stiffness 
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The  relationships  between  the  stress/moment  resultants  and  the  mid- 
surface  strains  are 
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(31) 


(32) 


(33) 


he  equations  of  motion  for  the  equivalent  continuum  plate  model  are 


The  equati 
[using  Eqs.  (20)] 
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where  q^,  q2,  and  are  the  external  load  components  and  M ^ ,  Mj,  and  M  3 
are  the  external  moment  components  in  the  x^,  x2 ,  and  directions, 
respectively. 
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p  .  dx„ 

'U  3 


PQ  x3  dx3 


Pq  x2  ^*3 


An  example  poblem  is  discussed  in  Section  IV. 


IV. 


EXAMPLE  PROBLEM 


A.  EQUIVALENT  STIFFNESSES 

In  order  co  assess  the  accuracy  of  using  the  equivalent  continuum  -ode'., 

the  methodology  was  applied  to  a  hypothetical  space  platform  (Fig.  3)  composed 
of  an  assembly  of  tetrahedral  trusses.  The  extensional,  coupling,  and  bending 
stiffnesses  are  determined  from  Eqs.  (28)  through  (30)  following  the  procedure 
given  by  Heki  (Ref.  29): 

Extensional  Stiffness 
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^Heki,  K.,  "On  the  Effective  Rigidities  of  Lattice  Plates,"  Recent 
Researches  of  Structural  Mechanics,  Tokyo,  1968. 
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The  density  parameters  for  the  equivalent  continuum  model  are  obtained 
from  Eqs.  (19)  and  (35): 
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With  the  help  of  Eqs.  (31) through  (33),  and  Eqs.  (36)  through  (38),  the 
equations  of  motion  (34)  are  rewritten  in  terms  of  the  displacements  and 
rotation  components  as  follows: 
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In  the  present  example  problem,  the  cross-sectional  areas  of  the  mem¬ 
bers  in  the  top  and  the  bottom  layers  are  assumed  to  be  identical  and  all  the 
members  lengths  are  equal,  i.e. 
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The  consequence  of  these  assumptions  is  that  p^ 


and  all  B.  's  are  zero. 
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Also  from  Eqs.  (38) 
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The  stiffness  and  mass  properties  of  the  equivalent  continuum  model  are 
marized  in  Table  4. 


Table  4.  Equivalent  Continuum  Model  Properties 


The  equations  of  motion  (40)  through  (45)  now  simplify  as 
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The  accuracy  of  the  equivalent  continuum  model  developed  is  illustrated 
by  the  following  free  vibration  problem  in  which  comparison  of  the  natural  fre¬ 
quencies  obtained  from  the  continuum  model  is  made  with  those  obtained  directly 
from  the  actual  truss  structure  solution.  For  the  free  vibration  case 

q.  s  0  i  -.1,2,3 

and 

M.  a  0 

i 
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Differentiating  Eqs.  (49)  and  (50)  with  respect  to  x^  and  x0 ,  respectively 
and  adding  them  together  vield 


/£*. 

a20  \ 

-  A 

9  . 

a“u3\ 

Ui 

•4) 

*3131 

\ »«; 

ax2/ 

A31 3 1^ 


7>  £jk 
at' 


(52) 


where 


d<P,  d<P0 
<P  H  ~  +  ~ 

axj^  dx2 


Equation  (48)  can  be  rewritten  as 

,,2  °  ,2  0 

/  a  113  a  U3  \ 

A3 131  (  2  +  2  )  +  A313l  ^  =  Pi 

\  0*2  / 


2  0 
a  U3 

0  ~2 
at 


or 


<P 


-  2  0 
P0  au3 


where  operator 


A3131  at2 


Jl  +  Jl 

ax2  ax2 


2  0 

V  u. 


Substitution  of  *l>  from  Eq.  (53)  gives 


(53) 


(54) 


P0  *2 


A3 131  at2 


2  0 

0  -  a  u3 
U3  P0  - 2 

at 


=  0 


(55) 


Equation  (55)  includes  the  rotary  inertia  as  well  as  transverse  shear  deforma¬ 
tion  terms. 
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If  Che  rotary  inertia  terms  are  omitted  from  Eq .  (55),  it  reduces  to 


*0  d2  \ 

A3 131  at2  / 


2  0 
T  u3 


d2  A 

^0  2 

at 


o 


(56) 


If  the  transverse  shear  deformation  is  neglected,  but  the  rotary  inertia  terms 
are  retained,  Eq.  (55)  reduces  to 


D  V 


-  a2  \ 

P2  ,2 
at 


) 


2  0 
V  u3 


2  0 

a  u3 

77 


(57) 


Finally,  if  both  the  transverse  shear  deformation  and  rotary  inertia  terms  are 
omitted,  Eq.  (55)  reduces  to  the  classical  plate  equation 


„  JO 
DV  u3 


(58) 


The  solution  of  Eqs.  (55),  (56),  (57),  or  (58)  for  appropriate  boundary  condi¬ 
tions  shall  yield  the  natural  frequencies  for  the  equivalent  continuum  model. 

The  significance  of  including  the  rotary  inertia  and  trans/erse  shear 
terms  has  been  discussed  by  Reissner  (Ref.  30),  and  Mindlin  (Refs.  31,  32). 


^Reissner,  E.f  "The  Effect  of  Transverse  Shear  Deformation  on  the  Bending 
of  Elastic  Plates,"  Journal  of  Applied  Mechanics,  Vol.  12,  June  1945. 

■^Mindlin,  R.  D. ,  "Influence  of  Rotary  Inertia  and  Shear  on  Flexural  Motions 
of  Isotropic,  Elastic  Plates,"  Journal  of  Applied  Mechanics,  Vol.  18,  No.  1, 
March  1951. 

■^Mindlin,  R.  D. ,  A.  Schacknow,  and  H.  Deresiewicz,  "Flexural  Vibrations  of 
Rectangular  Plates,"  Journal  of  Applied  Mechanics,  Vol.  23,  No.  3,  September 
1956. 
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Ic  is  known  that  both  effects  serve  to  decrease  the  computer  frequencies 

present  study,  the  natural  frequencies  are  computed  from.  Eqs.  (55;  through 
(58)  and  compared  with  those  obtained  from  solving  the  original  truss 
structure . 

As  a  first  example  case,  simply  supported  boundary  conditions  are 
assumed.  It  is  possible  to  obtain  the  continuum  model  closed  form  solution 
for  such  boundary  conditions.  These  boundary  conditions  may  be  stated  (see 


u 

Assuming  displacement  as 


A  s  in 


2  iuc 
—  e 


where  A  is  a  constant,  m  and  n  are  integers,  and  u>  is  circular  frequency.  It 
is  clear  chat  the  boundary  conditions  are  automatically  satisfied  by  the 
displacement  function.  The  frequency  U)  may  now  be  determined  by  substituting 
u°  into  Eqs.  (55)  through  (58). 

Case  1.  Rotary  inertia  and  transverse  shear  terms  included.  The 
characteristic  equation  in  this  case,  obtained  from  Eq.  (55)  is 


P  P 

y0r2  4 
: -  u 


Case  2.  Transverse  shear  term  included,  rotary  inertia  term  omitted. 
The  natural  frequencies  in  this  case,  obtained  from  Eq.  (56),  are  given  by 


Case  3.  Rotary  inertia  term  included,  transverse  shear  term  omitted. 
Ip  this  case,  the  natural  frequencies,  obtained  from  Eq.  (57),  are  given  by 
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Case  4.  Classical  plate  solution.  The  natural  frequencies  in  this 
case,  obtained  from  Eq .  f  53 ) ,  are  given  by 


B.  NUMERICAL  RESULTS 


The  following  geometric  and  material  properties  were  assumed  for  the 
tetrahedral  truss: 

Case  A 


A  *  A  =  50  x  10  ^  m2 

a 

£  -  6  m 

P  »  2  7  68  kg/m3 

E  *  71.7  x  106  kPa 

Case  B 

A  =  50  x  10'6  m2 

A  *  2 . 5  x  10  ^  m2 

a 

£  =  6  m 

P  -  2768  kg/m3 

71.7  x  106  kPa 


E  = 


The  natural  frequencies  for  the  tetrahedral  truss  with  6,  8,  12,  and  13 
repeating  elements  in  each  direction  were  computed  using  the  NASTRAN  computer 
program.  Equations  (60)  through  (63)  are  used  for  computing  the  natural  fre¬ 
quencies  of  the  equivalent  continuum  model.  If  N  is  the  number  of  the  repeat¬ 
ing  tetrahedral  truss  elements  in  each  direction,  plate  dimensions  a  and  b  are 
given  by 


a  =  Ni 


b 


N  f 


The  results  are  summarized  in  Tables  5  through  12  and  in  Figs.  8  through  15. 
The  finite  element  solution  of  the  truss  is  represented  by  FEM  in  the  tables, 
whereas  I  through  IV  represent  the  continuum  solutions  with  varying  degrees  of 
approximations.  These  are  defined  as  follows: 

I  =  Transverse  shear  and  rotary  inertia  terms  included 
II  =  Transverse  shear  terms  only  included 
III  *  Rotary  inertia  terms  only  included 
IV  =  Classical  plate  solution 


The  undeformed  shape  of  a  truss  with  12  repeating  tetrahedral  elements 
(orthographic  view)  is  shown  in  Fig.  16,  whereas  its  projections  on  the  X-Z 
plane  and  the  Y-Z  plane  are  shown  in  Figs.  17  and  18,  respectively.  The  cor¬ 
responding  views  for  the  first  four  mode  shapes  of  the  same  truss  structure 
are  shown  in  Figs.  19  through  30. 
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Table  5.  Comparison  of  Natural  Frequencies 


First  Mode 


CASc.  A 


- r 

i 

6  ! 

1  1 

8 

■  - ■—  —  - -  — 

12 

— 

18 

13. 67 

8.86 

1 

4.51 

2.18 

13.26 

8.28 

4.03 

1.88 

13.44 

8.37 

4.06 

1.88 

16.94 

9.  69 

4.36 

1.95 

17. 63 

9.92 

4.41 

1.96 

hi 


Number  of  Bays,  N 


Fig.  8.  Comparison  of  Natural  Frequencies  -  First  Mode 


Table  6.  Comparison  of  Natural  Frequencies  -  Second  Mode 


CASE  A:  A .  =  A 

a 


Table  8.  Comparison  of  Natural  Frequencies  -  Fourth  Mode 


CASE 

A:  A  , 

a 

A 

6 

8  | 

1 

1 

12 

1 

31.80 

23.  64 

12.86 

35.40 

24.16 

13.25 

35.80 

24.51 

13.44 

61.10 

36.41 

16.94 

70.  51 

39.66 

17.63 

Table  9.  Comparison  of  Natural 


59 


of  Bays,  N 


Frequencies  -  First  Mode 


Table  11.  Comparison  or  Natural  Frequencies  -  Third  Mode 


CASE  Sr  A  *  0.05. 

a 


FEM 

8.92 

6.  68 

4.36 

2.75 

I 

9.14 

6.79 

4.41 

2.78 

11 

9.14 

6.79 

4.41 

2.78 

III 

51.07 

30.23 

13.98 

6.33 

IV 

57.88 

32.  56 

14.47 

6.43 

Number  of  Bays,  N 

Fig.  14.  Comparison  of  Natural  Frequencies  -  Third  Mode 


Table  12.  Comparison  of  Natural  Frequencies  -  Fourth  Mode 


CASE  B:  A  *  '3.0  5A 

d 


FEM 

10. 63 

8.04 

5.33 

I 

11.15 

8.30 

5.43 

3.48  i 

i 

II 

11.15 

8.30 

5.43 

3.48 

III 

71.60 

43.16 

20.29 

9.26 

IV 

8  6.30 

47.98 

21.33 

9.48 

\ 

\ 


Number  of  Bays,  N 


Fig.  15.  Comparison  of  Natural  Frequencies  -  Fourth  Mode 


Vibration  Modes 

Tetrahedral  Truss  Structural  Model 
ly  Supported  Boundary  Conditions 


Vibration  Modes 


Undeformed  Shape  -  X-Z  Plane  Projection 
(Simply  Supported  Boundary  Conditions) 


Fig.  18.  Undef ort.ied  Shape  -  Y -7.  Plane  Projection 
(Simply  Supported  Boundary  Conditions) 


Vibration  Modes 

12  x  12  Tetrahedral  Truss  Structural  Model 
Simply  Supported  Boundary  Conditions 


First  Modt  -  Ortho^raphi  c  View 
(Simoly  Supported  Boundary  (.omniums 


Vibration  Modes 

12  x  12  Tetrahedral  Truss  Structural  Model 


N 
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■*• A-Ih 


Fig.  20.  First  Mode  -  X-Z  Plane  Projection 

(Simply  Supported  Boundary  Conditions) 


(Simply 


s*  E 


x  ^ 

CXI 


I 


Fig.  23.  Second  Mode  -  X-7,  Plane  Projection 

(Simply  Supported  Boundary  Conditions) 


Vibration  Modes 

12  x  12  Tetrahedral  Truss  Structural  Model 


Fig.  24 .  Second  Mode  -  Y-Z  Plane  Projection 

(Simply  Supported  Boundary  Conditions) 


Vibration  Modes 

12  x  12  Tetrahedral  Truss  Structural  Model 
Simply  Supported  Boundary  Conditions 


Vibration  Modes 

12  x  12  Tetrahedral  Truss  Structural  Model 
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Fig.  26.  Third  Mode  -  X-Z  Plane  Projection 

(Simply  Supported  Boundary  Conditions) 


Pig-  27.  Third  Mode  -  Y-Z  Plane  Projection 

(Simply  Supported  Boundary  Conditions) 


Vibration  Modes 

12  x  12  Tetrahedral  Truss  Structural  Model 
Simply  Supported  Boundary  Conditions 
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Fig.  29.  Fourth  Mode  -  X-Z  Plane  Projection 

(Simply  Supported  Boundary  Conditions) 


a> 

■§  c 

?  S 

TO  -O 

i=  c 

3  o 

o  <-> 

Si  £■ 

■§«^l 

5  w  c 
«/i  3 
r-  3  O 
O  >  00 

-.3  !i  -O 

£  TO  .22 

_o  »~  J- 

S  *p  o 
>  ®  a 
-C  a- 

TO  3 

*“  C/O 
-*-< 

a>  ^ 

t—  -S' 

CM 

—  t 

X  In 


Fig.  30.  Fourth  Mode  -  Y-Z  Plane  Projection 

(Simply  Supported  Boundary  Conditions) 


Numerical  results,  presented  for  the  simply  supported  tetrahedral  truss,  ;  : 

I  : 

demonstrate  the  accuracy  obtainable  from  the  equivalent  continuum  model.  It  is 
observed  that,  in  most  instances,  the  error  is  less  than  3  percent.  Further¬ 
more,  the  four  solutions,  I  through  IV,  demonstrate  the  significance  of  includ¬ 
ing  the  transverse  shear  and  rotary  inertia  terms  in  formulating  the  equivalent 
plate  problem.  It  is  observed  that  Solution  I,  which  includes  both  the  trans-  1 

verse  shear  and  rotary  inertia  effects,  is  very  good  even  when  the  number  of  ; 

basic  repeating  truss  elements  is  small.  The  effect  of  transverse  shear,  which 
is  quite  significant  for  a  truss  with  a  small  number  of  repeating  elements,  . 

becomes  less  pronounced  as  the  number  of  these  elements  increases.  This  is 
shown  in  Figs.  8  through  15.  The  rotary  inertia  effect,  on  the  other  hand,  is 
important  only  for  higher  frequency  responses  and  for  a  smaller  number  of 
repeating  elements.  Yet,  its  effect  is  much  less  pronounced  than  that  of  the 
transverse  shear  (see  Tables  5  through  12).  Furthermore,  i:  is  also  observed 
that  the  effect  of  transverse  shear  is  much  more  significant  for  the  case  in 

which  A,  =>  0.05A.  Finally,  the  mode  shapes  shown  in  Figs.  16  through  30  exhi- 
d 

bit  striking  similarities  to  those  of  a  plate  with  simply  supported  boundary 
conditions. 

4 

In  the  second  example  case,  free-free  boundary  conditions  are  assumed. 

The  natural  frequencies  for  the  truss  are  determined,  as  in  the  first  case, 
using  the  NASTRAN  computer  program.  Since,  for  the  free-free  boundary  condi¬ 
tion,  the  closed  form  solution  is  not  readily  available  for  the  equivalent  con¬ 
tinuum  plate,  the  natural  frequencies  are  determined  using  the  NASTRAN  computer 
program.  A  procedure  given  by  MacNeal  (Ref.  33)  is  followed  for  this  purpose. 


I 

^MacNeal,  R.  h. ,  "A  Simple  Quadrilateral  Shell  Element,"  Computer  and 
Structures ,  Vol.  8,  1978. 
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It  may  be  noted  here  that  the  finite  element  solution  of  a  free-free 
plate  accounts  for  transverse  shear  only.  The  natural  frequencies  computed 
for  the  truss  as  well  as  for  the  equivalent  continuum  model  are  given  in 
Tables  13  and  14. 


Table  13.  Comparison  of  Natural  Frequencies, 
Free-Free  Boundary  Conditions 


CASE  A: 


Frequency 


First  Mode 
*11 


Second  Mode 
f  21 


Third  Mode 
f12 


Fourth  Mode 
f22 


N 

6 

12 

Truss 

8.63 

2.56 

Continuum 

Model 

9.29 

2.66 

Truss 

12.47 

3.72 

Continuum 

Model 

13.81 

3.87 

Truss 

15.43 

4.91 

Continuum 

Model 

18.56 

5.50 

Truss 

18.80 

6.11 

Continuum 

Model 


20.07 


6.35 


Table  14.  Comparison  of  Natural  Frequencies 
Free-Free  3oundarv  Conditions 


CASE  B: 


F  requency 


First  Mode 

Truss 

*11 

Continuum 

Model 

Second  Mode 

Truss 

*21 

Cont inuum 
Model 

Third  Mode 

Truss 

*12 

Continuum 

Model 

Fourth  Mode 

Truss 

*22 

Continuum 

Model 

The  associated  mode  shapes  (orthographic  view,  projections  on  the  X-Z  and  Y-Z 
planes)  are  shown  in  Figs.  31  through  42.  A  good  comparison  is  once  again 
obtained  between  the  truss  model  and  its  equivalent  continuum  model.  The  mode 
shapes  (Figs.  31  through  42)  are  also  similar  to  those  of  a  free-free  plate. 


Vibration  Modes 

12  x  12  Tetrahedral  Truss  Structural  Model 
Free-Free  Boundary  Conditions 


First  Mode  -  Orthographic  View 
(Free-Free  Boundary  Conditions) 


Vibration  Modes 

12  x  12  Tetrahedral  Truss  Structural  Model 
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Fig.  32,  First  Mode  -  x-Z  Plane  Projection 
(Free-Free  Boundary  Conditions) 
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Second  Mode  -  X-7.  Plane  Pr 
(Free-Free  Boundary  Condit 


Boundary  Conditions 


Vibration  Modes 

2  Tetrahedral  Truss  Structural  Model 


ird  Mode  -  Orthographic  View 
ree-Free  Boundary  Conditions) 


Vibration  Modes 

Tetrahedral  Truss  Structural  Model 
'ee-Free  Boundary  Conditions 


Boundary  Condit 


Vibration  Modes 


40.  Fourth  Mode  -  Orthographic  View 
(Free-Free  Boundary  Conditions) 
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Fourth  Mode  -  X-Z  Plane  Projection 
(Free-Free  Boundary  Conditions 


Fig.  42.  Fourth  Mode  -  Y-Z  Plane  Projection 
(Free-Free  Boundary  Qjnditions) 


V.  CONCLUSIONS 


A  methodology  is  presented  for  modeling  large  truss-type  structures 
based  on  the  concept  of  equivalent  continuum-  The  procedure  for  obtaining 
equivalent  effective  elastic  -nd  dynamic  properties  in  terms  of  the  material 
and  geometric  properties  of  the  truss  is  simple  and  straightforward.  A  general 
three-dimensional  equivalent  continuum  model  is  reduced,  in  the  present  study, 
to  a  two-dimensional  plate  through  simple  kinematic  assumptions.  The  platform 
and  the  boundary  conditions  of  the  equivalent  continuum  plate  model  simulate 
those  of  the  original  truss  structures.  The  effects  of  transverse  shearing 
strain,  rotary  inertia,  bending-extensional  and  inertia  coupling  are  included 
in  the  continuum  model. 

Numerical  results,  presented  in  this  study,  are  for  a  tetrahedral  truss 
with  simply  supported  and  free-free  boundary  conditions.  Whereas  a  closed  form 
solution  was  obtained  for  the  equivalent  continuum  plate  with  simply  supported 
boundary  conditions,  the  same  was  not  possible  in  the  case  of  a  plate  with 
free-free  boundary  conditions.  The  solution  in  the  latter  case  was  obtained 
using  the  NASTRAN  computer  program.  The  results  adequately  demonstrate  the 
accuracy  obtainable  from  the  continuum  model.  In  most  instances,  the  error 
is  less  than  3  percent.  Furthermore,  the  four  solutions,  I  through  IV,  demon¬ 
strate  the  significance  of  including  the  transverse  shear  and  rotary  inertia 
terms.  It  is  observed  that  the  plate  solution  I,  which  includes  both  trans¬ 
verse  shear. and  rotary  inertia  effects,  is  very  good,  oven  when  the  number  of 
basic  repeating  truss  elements  is  small. 

For  the  truss  with  a  small  number  of  repeating  elements,  the  effect  of 
transverse  shear  is  significant.  As  the  number  of  the  repeating  elements  in¬ 
creases,  the  transverse  shear  effect  becomes  less  pronounced  as  shown  in  Figs. 

8  through  15.  The  rotary  inertia  effect,  on  the  other  hand,  is  important  only 
for  higher  frequency  responses  and  for  a  smaller  number  of  repeating  elements. 
Yet,  it  is  much  less  pronounced  than  that  of  the  transverse  shear  (see  Tables  5 
through  12).  Furthermore,  the  effect  of  transverse  shear  becomes  much  more 
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significant  for  the  case  in  which  =  0.05A  as  compared  to  the  case  in  which 

A  .=  A.  The  mode  shapes  of  the  truss  with  simplv  supported  and  free-free 
d 

boundary  conditions  (see  Figs.  16  through  24)  exhibit  striking  similarity  with 
those  of  a  plate  with  similar  boundary  conditions.  The  main  conclusions  of  the 
study  are  summarized  as  follows: 

1.  The  technique  of  replacing  a  large  truss  structure  with  an  equivalent 
continuum  is  extremely  promising. 

2.  Good  agreement  between  the  discrete  and  its  equivalent  continuum  model 
solutions  can  be  obtained. 

3.  The  classical  plate  solution  may  only  be  suitable  for  simple  geometries 
with  large  numbers  of  repeating  modules.  For  relatively  small  numbers 
of  repeating  truss  elements,  the  transverse  shear  and  rotary  inertia 
effects  must  be  included  in  the  so-lution  of  the  equivalent  continuum 
model . 
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